Experimental quantum computing without entanglement 
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Entanglement is widely believed to lie at the heart of the advantages offered by a quantum computer. 
This belief is supported by the discovery that a noiseless (pure) state quantum computer must 
generate a large amount of entanglement in order to offer any speed up over a classical computer. 
However, deterministic quantum computation with one pure qubit (DQCl), which employs noisy 
(mixed) states, is an efhcient model that generates at most a marginal amount of entanglement. 
Although this model cannot implement any arbitrary algorithm it can efficiently solve a range of 
problems of significant importance to the scientific community. Here we experimentally implement 
a first-order case of a key DQCl algorithm and explicitly characterise the non-classical correlations 
generated. Our results show that while there is no entanglement the algorithm does give rise to other 
non-classical correlations, which we quantify using the quantum discord — a stronger measure of non- 
classical correlations that includes entanglement as a subset. Our results suggest that discord could 
replace entanglement as a necessary resource for a quantum computational speed-up. Furthermore, 
DQCl is far less resource intensive than universal quantum computing and our implementation in 
a scalable architecture highlights the model as a practical short-term goal. 



In contrast to the highly pure multi-qubit states re- 
quired for the conventional models of quantum comput- 
ing [U [5], DQCl employs only a single qubit in a pure 
state, alongside a register of qubits in the fully mixed 
state [3]. While this model is strictly less powerful than 
a universal quantum computer (where one can implement 
any arbitrary algorithm) it can still efficiently solve im- 
portant problems that are thought to be classically in- 
tractable. The application originally identified was the 
efficient simulation of some quantum systems [3] . Since 
then exponential speed-ups have been identified in esti- 
mating: the average fidelity decay under quantum maps 
[3] ; quadratically signed weight enumerators [5^ ; and the 
Jones Polynomial in knot theory j6j . Recently it has been 
shown that DQCl also affords efficient parameter estima- 
tion at the quantum metrology limit f7^. Furthermore, 
attempts to find an efficient way of classically simulating 
DQCl have failed [5]. These results provide strong evi- 
dence that a device capable of implenting scalable DQCl 
algorithms would be an extremely useful tool. 

Besides the practical applications, DQCl is also fasci- 
nating from a fundamental perspective. For example, 
it is straightforward to show that a model employing 
only fully mixed qubits offers no advantage over a clas- 
sical computer. It is therefore surprising that the addi- 
tion of only a single pure qubit offers such a dramatic 
increase in computational power. Furthermore, an im- 
portant quantum information result is that a pure state 
quantum computer can only offer an advantage over a 
classical approach if it generates an amount of entangle- 
ment that grows with the size of the problem being tack- 
led IHITD]. This supports the commonly held belief that 
entanglement lies at the heart of the quantum advantage. 
However entanglement is at most marginally present in 
DQCl [11]. The existence of efficient DQCl algorithms 
therefore provides strong evidence that a large amount of 
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FIG. 1: DQCl normalised trace estimation algorithm. Re- 
peated running of the circuit and measurement of qubit c 
yields estimates of the standard Pauli {X) and (Y) expecta- 
tion values, from which the normalised trace of the unitary 
operator Un can be derived (Tr[l7n]/2"). /„ is the n qubit 
identity. 

entanglement is not necessary to achieve a speed-up. In- 
stead it has been proposed that a requirement, common 
to both pure and mixed state models, is the generation 
of quantum correlations that are not fully captured by 
entanglement [8j [12]. It is now important to pursue a 
deeper understanding of these correlations. 

Following recent developments in experimental quan- 
tum computing [13j we present an implementation of 
DQCl in a linear optic architecture. We build on a 
previous demonstration in liquid state NMR |14j in two 
important ways: we are able to characterise the quan- 
tum correlations generated by our implementation; there 
are several known paths to efficient scalable linear optic 
quantum computing [S] 115} 116] and there is much recent 
and ongoing progress towards developing the necessary 
technology [TT] [THj [IH] . 

We present a first-order implementation of the DQCl 
algorithm for estimating the normalised trace of a unitary 
matrix [3] [51 [TTJ [T5] . This algorithm achieves an expo- 
nential speedup over the best known classical approach 
i.e. it requires exponentially fewer resources as the size of 
the unitary increases. It is thought highly unlikely that 
an efficient, but as yet unknown, classical approach can 
exist [11]. That DQCl can perform this task efficiently 
underpins the model's ability to solve the aforementioned 



range of important practical problems. 

Figure 1 shows the normalised trace estimation al- 
gorithm. The required input state is separable and 
consists of a single pure qubit c in the logical state 
|0)(0|=i{/i + Z}, and a register of n qubits in the com- 
pletely mixed state In/'^"- The circuit consists of the 
standard Hadamard gate jT] and a unitary ([/«) on the 
register, controlled by qubit c. The state of all qubits 
at the output of the circuit is: 
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where N—2" and /„ is the n qubit identity. The reduced 
state of qubit c (after performing a partial trace over the 
register) is given by: 
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Thus the normalised trace of J7„ is encoded in the co- 
herences of qubit c. This information can be retrieved 
by measuring the expectation values of the standard 
Pauh operators X and Y, since {X)= Re[Tr{U„) /2N] and 
{Y)=~lm[Tr{Un)/2N]. 

An expectation value is estimated by repeatedly run- 
ning the circuit. One can achieve a fixed accuracy e in 
this estimate with a number of runs L ^ ln(l/Pe)/e^, 
where Pe is the probability that the estimate is farther 
from the true value than e (Ref. [H])- That the ac- 
curacy does not scale with the size of the unitary and 
scales logarithmically with the error probability means 
that this is an efficient algorithm for estimating the nor- 
malised trace. In contrast, classical approaches suffer an 
exponential increase in required number of resources with 
the size of the unitary [TT] . 

We implement the first-order {n=l) case for the uni- 
tary: 
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In this case {X)={1+ cos9)/2 and (r) = (sin 6')/2. 

Our implementation is performed in an all optical ar- 
chitecture shown in Fig. 2 (see Methods). We encode 
quantum information in multiple degrees of freedom of 
single photons. Single qubit gates are realised determinis- 
tically using birefringent wave-plates. The required two- 
qubit gate is realised non-deterministically and measure- 
ment of single photons in the two output modes signals 
a successful run [T3] . 

The trace estimation algorithm is implemented over 
the range — 7r<6'<7r (Eqn. 3). Fig. 3a compares the ex- 
perimentally observed results with the ideal. The ideal 
case is calculated assuming perfect circuit operation and 
measured input states. We observe a high degree of cor- 
relation quantified by a reduced of 0.7 (real curve) and 
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FIG. 2: Experimental schematic. The controUed-Zg gate is 
implemented using a recently developed technique |13j . At 
the input and output qubits are encoded in the polarisation 
of single photons (|0)=|h), |1) = |v), Horizontal and Vertical). 
Coincident measurement of single photons at detectors Dl- 
2 signals a successful run of the gate. The interferometers 
in each mode allow the preparation of noisy (mixed) states, 
by introducing a path difference between the two arms. See 
Methods. 

1.2 (imaginary curve). Deviations from the ideal are due 
to imperfect circuit operation caused by: optical beam 
steering as 9 is varied; interferometric instability; and 
non-classical interference instability. There are several 
known paths to reduce these errors including moving to 
micro-optic systems [15]. 

Interestingly, the exponential speed-up offered by this 
algorithm is not affected by reducing the purity of qubit 
c llj. Consider replacing the initial state of this qubit 
with the mixed state ^{/i -I- ctZ}, where a now reflects 
the purity (p=[l + a'^]/2, 0<q;<1). At the output of the 
circuit the state of this qubit is now given by: 

_ r 0.5 aTr{Un)yN' 
~ aTr{Un)/N 0.5 



(4) 



The effect of mixture in qubit c is to reduce the expecta- 
tion values {X) and (Y) by a (Eqn. 2), thereby making 
it harder to estimate the normalised trace. In order to 
achieve the same fixed accuracy as before requires an in- 
creased number of runs L'^L/a^ . While this clearly adds 
an additional overhead, as long as a is non-zero, the al- 
gorithm still provides an efficient evaluation of the trace. 
Even access to the tiniest fraction of a single pure qubit 
is sufhcient to achieve an exponential speedup over the 
best known classical approach. 

Fig. 3b compares the experimentally observed algo- 
rithm results with the ideal, for Q!=0.58±0.02. We ob- 
serve a high correlation with the ideal, quantified by a 
reduced of 1.8 (real curve) and 2.0 (imaginary curve). 
The increased in this case (compared to the results in 
Fig. 3a) is due to a less favorable optical alignment and 
is not an intrinsic error associated with initialising c into 
a mixed state. The additional resource overhead in this 
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FIG. 3: Algorithm output. Real and imaginary parts of the 
renormalised trace measured for two values of a, over a range 
of d, Eqn. 4. Theoretical predictions are calculated using 
measured input states and assuming perfect circuit operation. 
We observe a of a, 0.95 and b, 2.0. (X) is estimated by 
counting the number of coincident photon pairs (A'') when 
projecting qubit c into the states |±)=(|0)±|l))/\/2, over 10 
sees. Then {X)=(iV+-iV_)/(A^++Af_). The same technique 
is used to estimate (y), but in this case we project into the 
states |±i) = (|0)±ill))/\/2. 

case is apparent in the amplitude reduction by a factor 
of OL compared with the results shown in Fig. 3a. 

Finally we analyse the correlations generated by the al- 
gorithm by performing full tomography of the two-qubit 
output state before the final measurement stage (Eqn. 1), 
using 36 (over-complete) measurement bases. This allows 
a reconstruction of the density matrix. Figure 4 shows 
two measures of non-classical correlations, derived from 
the experimentally reconstructed density matrices — the 
well-known tangle |20j and the lesser-known quantum 
discord [12 |^ The tangle is a complete measure 
of entanglement in two-qubit states and represents per- 
haps the most striking divergence from classical behav- 
ior. However, entanglement is not the only kind of non- 
classical correlation. A far stronger measure, that encom- 
passes entanglement and more, is given by the quantum 
discord. 

The discord is concerned with a fundamental charac- 
teristic of classical systems — that their information con- 
tent is locally accessible and can be obtained without 
perturbing the state for an independent observer [55]. 
If the discord is zero there exists a local measurement 
protocol under which all the state information can be re- 
vealed without perturbing the state for observers who do 
not have access to the measurement results. If the dis- 
cord is non-zero then no such protocol exists. For pure 
states, discord becomes a measure of entanglement — no 
other non-classical correlations can be distinguished in 



this case. However, for mixed states the discord captures 
more non-classical correlations than entanglement |12j . 

The results show that, to within experimental error, 
our implementation does not give rise to any entangle- 
ment. However, in general it does generate quantum dis- 
cord. We observe a high degree of correlation between 
the theoretical and measured discord values, quantified 
by a of 1.6. These results are consistent with recent 
theoretical work [12] which predicts that, although the 
entanglement is generally zero for arbitrary instances of 
this algorithm, discord is consistently present. 

For our implementation the discord is zero in two dis- 
tinct cases 0={O, ±7r}, corresponding to our controUed-Zg 
gate implementing the identity and a controlled-Z7r , re- 
spectively. Both of these gates, alongside the Hadamard, 
are members of the Clifford group [T]. In these cases 
the entire state evolution is implemented by gates from 
the Clifford group. Furthermore, the algorithm involves 
preparing the input in a mixture of logical basis states 
and measurement of observables in the Pauli group [1]. 
Under these conditions the Gottesman-Knill theorem 
tells us that the entire algorithm can always be efficiently 
simulated on a classical computer [Ij (see supplementary 
material). In contrast, for all other values of Q the ac- 
tion of the controUed-Ze gate is responsible for a non- 
Clifford group evolution. There is no known classical 
method to efficiently simulate an arbitrary size algorithm 
that evolves in this way — thereby allowing for a quantum 
speed-up. Furthermore, it is straightforward to show that 
any order implementation of the normalised trace esti- 
mation algorithm that is composed entirely of gates from 
the Clifford group produces a state with zero discord (see 
supplementary material) . These results provide evidence 
of a link between discord and the potential for a DQCl 
speedup. An important path for further research is to 
determine whether all DQCl circuits that do not gen- 
erate discord can be efficiently simulated on a classical 
computer. Such a result would provide strong evidence 
that the discord is a more accurate measure of the re- 
sources required for a quantum speed up, than entangle- 
ment. Besides the fundamental interest, this could have 
implications in the many burgeoning quantum computing 
architectures where environmental decoherence presents 
a significant obstacle to universal pure state quantum 
computing. 

We have generated, and fully characterised, quantum 
correlations that lie beyond entanglement and provided 
evidence that they enable efficient quantum computa- 
tion. It is now of interest to explore quantum discord 
in the context of other active research areas, such as 
'non- locality without entanglement' [531 Hi]- while 
the two-qubit states of interest in these works are not 
entangled they all have non-zero discord, signifying the 
presence of quantum correlations. 

The authors wish to thank C. Caves, A. Datta, A, 
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FIG. 4: Non-classical correlations generated by the algo- 
rithm. Discord and tangle derived from experimentally re- 
constructed density matrices measured at the algorithm out- 
put for Q=0.997±0.003 (Fig. 3a). The discord is calculated 
by optimising over all 1-qubit local projective measurements 
on qubit c, Fig 1 (see supplementary material). Theoretical 
predictions are calculated using measured input states and 
assuming perfect circuit operation. The observed for the 
discord is 1.6. 
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METHODS 

Photons are generated via spontaneous parametric down 
conversion of a frequency-doubled mode-locked Ti: Sapphire 
laser (820nm^410nm, Ar=80fs at 82MHz) through a Type-I 
2mm BiBaOe crystal. Photons are filtered by blocked inter- 
ference filters at 820±1.5nm; collected into two single-mode 
optical fibers; injected into free-space modes c and r (Fig. 2); 
detected using fibre-coupled single photon counting modules 
(D1-D2). With lOOmW at 410nm we measure a two-fold co- 
incidence rate at the output of the optical circuit of approx- 
imately 100 s~^. The required CNOT is implemented using a 
standard technique involving non-classical interference at a 
partially polarising beamsplitter and projective measurement 
|251 1261 127j . State dependent loss is used to rebalance state 
amplitudes [251 1261 [27] . To improve count rates, we achieve 
correct balance by pre-biasing the c input state [251 1261 \T7\ . 
The algorithm success probability is 1/12. We realise the re- 
quired interferometers using calcite beam displacer pairs [28] 
(Fig. 2). 

Each qubit goes through a polarisation interferometer 
(Fig. 2). We control the level of mixture in a qubit by al- 
tering the path-length difference between the arms of the 
corresponding interferometer. A path difference that is 
greater than the photon coherence length results in a fully 
decohered — that is, a fully mixed — photonic qubit. By tun- 
ing the path difference between zero and the photon coherence 
length we can accurately control the level of mixture in the 
qubit between zero and maximum. We alter the path differ- 
ence by rotating one calcite beam displacer of a pair about 
an axis perpendicular to the plane defined by the two paths. 

All error bars are calculated through photon counting 



uncertainty described by poissonian statistics. We use the 
standard definition for a reduced calculation, allowing for 
three degrees of freedom (in our implementation both the real 
and imaginary parts of the trace are simple trigonometric 
functions defined by an amplitude, frequency and phase, 
Eqn. 3). 



SUPPLEMENTARY INFORMATION. 



The Gottesman-Knill theorem for mixed states. 

The Gottesman-Knill theorem applies to algorithm input 
states prepared in the logical basis. In the DQCl model 
the input is a mixture of logical basis states [l]. However, 
it is straightforward to show that such an input state can 
be prepared from a pure logical state using only a number 
of additional Clifford group gates and ancilla qubits that is 
linear in the size of the input state. As such the DQCl input 
state satisfies the conditions under which the theorem is valid. 

Definition of quantum discord for two-qubits. 

The definition of the mutual information for a bipartite den- 
sity matrix is [nil2|: 



(5) 



where H{p) is the well-known Von Neumann entropy of the 
state [T] p. An alternative definition is given by: 



J{r:c)^H{pr) - H{p,r\c), 



(6) 



where H{pcr\c) is the extension of the classical conditional en- 
tropy to the quantum case '17; . This is obtained by minimis- 
ing the average entropy of the subsystem r, over all possible 
projective conditional measurements on c: 



H{pcr\c) = min^pi_f/'(prin. 



where pi=Tx{Ilipcr^i), and p^p^ = Trc{J\-ipcr^i)/pi 
discord is defined as the difference [12[ : 

V{r,c) = l{r:c)-J{r:c). 
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Notice that J'(r-:c) is not symmetrical by inversion of c and r, 
therefore, in general, discord is directional: T>{r,c) 7^ T>{c,r)\ 
we might not be able to detect quantum correlation when 
conditioning on measurements of one partition, while they 
arise when considering the inverse case. It is straightforward 
to show that states admitting a diagonal representation in 
a local basis have bidirectionally vanishing discord — they 
contain only classical correlations. 

Proof that DQCl Clifford evolution generates 
no discord. 

The DQCl input state can be written in the form 
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which is clearly diagonal in the logic basis, hence it has zero 
discord in both directions. The action of Clifford group gates 
is to map Pauli matrices into Pauli matrices: if we indicate 
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the unitary action of the circuit by W, the input state pi„ is 
transformed into: 
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where ar'' refers to the i-th qubit, and r can take the values 
r = 0,1,2,3 labeling the Pauli matrices: ao — I, ai — 7i, 
02 = X, and aa = Y. The state ( 10 1 is locally equivalent to 
the state: 

P'=^(l«"+^ + 0-i")> (11) 

where the index s can take only the values s = 0, 1; thus p' 
admits a diagonalisation in a local basis. This state is purely 
classically correlated, hence its discord is zero. Consequently 
Pout, which is obtained from p' with local rotations, must 
have zero discord. 
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